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Abstract 



Various methods of searching for supersymmetric dark matter are sensitive 
to WIMPs with different properties. One consequence of this is that the 
phenomenology of dark matter detection can vary dramatically in different 
supersymmetric breaking scenarios. 

In this paper, we consider the sensitivities to supersymmetric dark matter 
of different detection methods and techniques in a wide variety of super- 
symmetric breaking scenarios. We discuss the ability of various astrophysical 
experiments, such as direct experiments, gamma-rays satellites, neutrino tele- 
scopes and positron and anti-proton cosmic ray experiments, to test various 
supersymmetry breaking scenarios. We also discuss what information can 
be revealed about supersymmetry breaking by combining results from com- 
plementary experiments. We place an emphasis on the differences between 
various experimental techniques. 
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I. INTRODUCTION 



An enormous body of evidence has accumulated in the favor of cold dark matter. This 
body of evidence includes observations of galactic clusters and large scale structure [1], 
supernovae [2] and the cosmic microwave background (CMB) anisotropics [3,4]. Recently, 
WMAP has provided the most detailed information on the CMB to date, quoting a total 
matter density of fl^h^ = O.lSSlHog 

[4]. Furthermore, data from WMAP and other prior 
experiments indicate a considerably smaller quantity of baryonic matter [4,5]. At the 2-a 
confidence level, the density of non-baryonic, and cold, dark matter is now known to be 
QcBMh^^O.nstlill [4]. 

A compelling candidate for the cold dark matter is naturally provided by supersymmetry 
[6]. In supersymmetric models which conserve R-parity [7], the lightest supersymmetric 
particle (LSP), is stable. Furthermore, in many supersymmetric models, the LSP is the 
lightest neutralino, a mixture of the superpartners of the photon, Z and neutral higgs bosons, 
and is electrically neutral, colorless and a viable dark matter candidate. If such an particle 
were in equilibrium with photons in the early universe, as the temperature decreased, a 
freeze-out would occur leaving a thermal relic density. The temperature at which this occurs, 
and the density which remains, depends on the annihilation cross section and mass of the 
lightest neutralino. It is natural for supersymmetry to provide a dark matter candidate with 
a present abundance similar to those favored by the WMAP experiment. 

Many methods have been proposed to search for evidence of neutralino dark matter. 
In addition to accelerator searches for supersymmetry [8], direct and indirect dark matter 
detection experiments have been performed. Direct dark matter searches attempt to measure 
the recoil of dark matter particles scattering elastically off of the detector material [9,10]. 
Indirect dark matter searches hope to observe the products of dark matter annihilation 
including gamma-rays [11-15], neutrinos [16], positrons [17-19] and anti-protons [20]. 

Supersymmetry with low scale soft breaking has been established as an excellent candi- 
date for physics beyond the standard model. At low energy, a usual practice in studying 
supersymmetry phenomenology is to parameterize the supersymmetry breaking by putting 
in explicit soft supersymmetry breaking terms in the effective Lagrangian. In principle, all 
of those terms should be derived from some supersymmetry breaking and mediation mech- 
anism which is embedded in the fundamental theory. Therefore, supersymmetry breaking 
mechanisms provide the crucial link between the fundamental theory and the low energy pa- 
rameters. Knowledge of supersymmetry breaking will provide us with important information 
about the physics at the high energy scale. We would like to extract as much information 
as possible about how supersymmetry is broken from our low energy experiments. Needless 
to say, the identity and properties of the LSP would carry invaluable information about the 
underlying supersymmetry breaking mechanism. In this paper, we attempt to evaluate the 
prospects for direct and indirect detection of neutralino dark matter in a variety of supersym- 
metry breaking scenarios. We find distinctive differences in the dark matter phenomenology 
for various breaking mechanisms. Information from cold dark matter experiments may be 
very useful in providing hints about which supersymmetry breaking scenario is manifest in 
nature. 

This paper is organized as following: First we briefiy review various direct and indirect 
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detection methods for cold dark matter. We then review the rclcvcnt information regarding 
the different supersymmetry breaking scenarios that we study in this paper. Finally, we 
present our results and comparison of different supersymmetry breaking scenarios. 

II. DIRECT DETECTION 

If the dark matter halo of our galaxy is made up of TeV-scale particles, then millions 
of such particles travel through each square meter each second in our galaxy. Despite this 
large rate, dark matter particles likely have very small cross sections, making observation of 
such particles very difficult. 

Proposed as early as twenty years ago [9] , numerous experiments have been developed in 
an effort to directly detect dark matter particles. In this paper, we will compare theoretical 
predictions to recent experimental results and to the reach of future efforts. 

The ability of a direct detection experiment to observe a WIMP depends on that particle's 
elastic scattering cross section with nucleons in the detector's target medium. This quantity 
depends strongly on the composition of the lightest neutralino (gaugino verses higgsino) as 
well as the masses of the squarks mediating scattering processes. 

Limits or sensitivities are most frequently shown as a elastic scattering cross section 
verses WIMP mass. Since the lightest neutralino is a Majorana fermion, it only has axial- 
vector and scalar interactions with nuclear matter. For heavy nuclei, typically A > 20, the 
scalar (or spin independent) interaction dominates. Therefore, we focus on the limits and 
reaches for scalar interaction cross sections. For each model, we assume the measured local 
dark matter density and velocity distribution [21]. 

There are a large number of direct detection experiments which have been developed or 
are planned. Currently, the strongest sensitivities have been attained by the CDMS [22], 
ZEPLIN-I [23] , and Edelweiss [24] experiments. Additionally, the DAMA collaboration has 
claimed an observation of an annual modulation in their rate consistent with a ~ 50 GeV 
WIMP with (T^i\f ~ 10~^ pb [25]. The DAMA result has been highly controversial and the 
entire region claimed by the experiment is now excluded by other experiments. 

The next round of direct experiments include CREST-II [26], CDMS (Soudan site) [27], 
Edelweiss II [28] and ZEPLIN-II, III, and IV [29]. These experiments will be capable of 
studying WIMPs with considerably lower cross sections in the relatively near future. Later, 
experiments such as XENON [30] or GENIUS [31] may probe even further. 

III. INDIRECT DETECTION 

Several different methods have been pursued to search for cold dark matter particles 
indirectly. In these approaches, experiments are designed to look for annihilation products 
from WIMP annihilations, such as neutrinos from the Sun or Earth, gamma-rays from the 
center of the galaxy, and positrons or anti-protons from the galactic halo. 
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A. Gamma-Rays From The Galactic Center 



In regions where the dark matter density is very large, the annihilation rate may become 
large enough to provide observable fluxes of high energy gamma-rays (and also neutrinos). 
The spectral [32] and angular [33] features of such a source have been studied for decades. 
The center of our galaxy may be such a region. Satellite-based experiments, such as EGRET 
and GLAST, are capable of observing gamma-rays at the energies suitable for this purpose. 

The rate at which neutrahnos annihilate near the galactic center depends strongly on 
the halo profile considered. At this time, the distribution of galactic dark matter, especially 
near the galactic center, is not well known. The halo models commonly discussed in this 
debate include distributions with low density cores, high density cusps or even higher density 
spikes. 

N-body simulations appear to predict models with cuspy distributions, such as the 
Navarro, Prenk and White (NFW) or Moore, et. al. profiles [34]. In such models, the 
dark matter density increases as, p oc 1/r'^, near the central region. 7 is 1.0 in the NFW 
profile and 1.5 in the Moore, et. al. profile. Arguments have been made, however, that 
observations of galactic rotation curves favor fiat density core profiles [35] , although others 
claim that these observations do not preclude the presence of a cusp at the center of our 
galaxy [36] . In the case of a core profile, no observable gamma-ray signal would be produced 
in dark matter annihilations near the galactic center region. 

Models with substantial central density spikes have been discussed a great deal recently 
[37] . It has been argued that density spikes are generated naturally as a result of adiabatic 
accretion of matter into the central galactic black hole. If density spikes form in the center 
of galaxies, it would be possible for dark matter annihilation rates to occur near the center 
of our galaxy which were considerably larger than those predicted for cuspy or other halo 
profiles. 

There are two signatures for the detection of gamma-rays from the galactic center which 
we will consider in this paper. The first is line emission produced from the processes x^x'^ — > 
77 and x'^x° ~^ 1^- The second is continuum emission from all gamma-ray producing 
processes, dominated by 7r° decay created in the fragmentation of quarks [12]. Although the 
continuum emission lacks the distinctive spectral features of the line emission, the number 
of photons produced is considerably greater. Since there is no direct coupling between 
the neutralino and the photon, all line producing processes must be mediated by some 
charged particle loop (such as SM fermions, charginos, W-bosons, etc.). Therefore, all hue 
emission processes are necessarily loop suppressed. For all of the diagrams contributing to 
line emission at the one loop level, see [13]. 77 and 7Z line emission will give rise to very 
distinct feature in the spectrum (since the end products will have a very narrow energy 
distribution). The position of these lines can also be used to measure the mass of the 
annihilating WIMP. This is probably the only way to accurately measure the WIMP's mass 
outside of a coUider. 
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B. Neutrinos From The Sun 



In addition to the galactic center region, high densities of dark matter may accumulate in 
the center of less distant objects such as the Sun and Earth. In these deep gravitational wells, 
the annihilation rate of WIMPs can be large, resulting in observable fluxes of annihilation 
products, such as neutrinos [38]. High-energy neutrino telescopes, such as AMANDA [39], 
ANTARES [40] or next generation IceCube [41], are designed, in part, for this task. For a 
review of high-energy neutrino astronomy, see Ref . [42] . 

The flux of neutrinos from neutralino annihilations in the Sun or Earth is a function 
of the annihilation cross section of the WIMP, the WIMP-nucleon elastic scattering cross 
section and the ratios of the various WIMP annihilation modes. For most SUSY models, the 
annihilation rate reaches equilibrium with the capture rate in the Sun, and the dependence on 
the annihilation cross section is removed. Typically the flux of neutrinos, and the associated 
event rate, from WIMP annihilation in the Earth is much far smaller than from the Sun, so 
we will consider only the rate from the Sun in this paper. 

Neutrinos can be produced in neutralino annihilation by several processes including 
^ tt, x°X° ^ bb, x°X° ^ cc, x°X° ^ ZZ, x°X° ^ W+W' and x°X° ^ r+r" 
[43]. Neutrinos are produced directly in the decays of r='='s, c and b quarks and gauge 
bosons, and indirectly through the decays of the b quarks and bosons created in top 
quark decays. The process(es) which dominates depends on the mass and composition of 
the lightest neutralino. In the fragmentation process, particles can lose substantial amounts 
of their annihilation energy before they decay and produce a neutrino. Furthermore, high- 
energy neutrinos can interact as they propagate through the Sun, thus degrading their flux 
[43,44]. 

C. Positrons and Anti-Protons From The Galactic Halo 

Annihilating halo neutralinos can also generate positrons and anti-protons. Unlike 
gamma-rays or neutrinos, charged cosmic rays do not point at their sources due to galactic 
magnetic flelds. By measuring the flux of positrons and anti-protons at Earth, it may be 
possible to observe signatures of annihilating neutralinos in nearby dark matter clumps. 

In 1994 and 1995, the High Energy Antimatter Telescope (HEAT) observed an excess 
of cosmic positrons, peaking around ~ 10 GeV [45]. This result was confirmed by another 
HEAT flight in 2000 [46] . Although it is uncertain whether this excess is the result of dark 
matter annihilation, it has been shown that halo positrons could be a signature of such a 
process [17-19]. 

Using a smooth halo profile, we calculate the positron fiux on Earth following the meth- 
ods of [19]. This includes not only the calculation of the injection rate of positions from 
neutralino annihilation, but also the propagation effects such as synchrotron and inverse 
Compton energy losses, which are very important. We find that for a smooth halo dis- 
tribution, the positron flux predicted is probably insufficient to explain the HEAT data. 
Any deviations from a smooth halo profile, however, enhance the positron rate, especially 
if clumps occur relatively nearby. We parameterize this effect with a single variable, the 
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positron boost factor, which is the ratio of the positron flux at Earth for a given halo profile 
to the positron flux for a smooth profile. 

The anti-proton flux can also be calculated [20]. Anti-protons can propagate much 
greater distances than positrons, however, often travehng up to ten or more kiloparsecs. This 
requires a separate boost factor parameter for anti-protons. The most relevant measurement 
of the anti-proton flux comes from the 1995 and 1997 BESS (Balloon Borne Experiment 
with Superconducting Spectrometer) data [47]. The anti-proton flux measured by BESS is 
1.27l[j:i X lO-^cm-^s-^sr-^GeV-^ in the range of 400 to 560 MeV. 



IV. SUPERSYMMETRY BREAKING 



The general framework to study the phenomenology of softly broken low energy super- 
symmetry is the Minimal Supersymmetric Standard Model (MSSM). The matter spectrum 
of the MSSM can be obtained basically by assigning a superpartner to each of the Standard 
Model fields. The only exception is that in the MSSM, we require two Higgs doublets. 
The MSSM superpotential is given by 

W = ea,[-KQ\Y:m^^ + H^.QtYl'D^j + HfLlY^^E'^, - ^iH^.H':], (1) 

in which Cab — —^ba and eu — 1, and the superfields are defined in the standard way 
(suppressing gauge indices): 



Q^ = 




U! = 




A— 






{El,, Eli 


Et = 




Hu = 


{Hu. Hu) 


Hd = 


{Hd, Hd), 



(2) 

with i,j — 1 ... 3 labeling family indices. The soft breaking Lagrangian, Csoft, takes the 
form (dropping "helicity" indices): 



1 



■'soft 



M^gg + M2W''W^ + MiBB 



+ eab[-hHlHl - H:Q\A^^.U] + H^dQ'iAd.fi'i + HfL\A,^.E'^^ + h.c] 
+ ml,\Hd\'' + mlJ\Hj\^ + g>|. .Qf 

+ ^>i./r + Urmlip'', + DTm%p'', + Erm%^W^. (3) 

The m^'s, usually called soft masses squared, are in general 3x3 Hermitian matrices. A are 
usually called trilinears which are 3x3 general complex matrices. In many supersymmetry 
breaking scenarios where A is universal, such as mSUGRA (which we will review later), it is 
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common to write A — Y ■ A. Mj's are the mass terms for the superpartners of the Standard 
model gauge bosons, gauginos. 

Notice that the MSSM by itself is not a model of supersymmetry breaking. It is a param- 
eterization of supersymmetry breaking by including general soft breaking terms explicitly. 
As a result, the MSSM has a large parameter space containing 124 parameters. Most of these 
are soft SUSY breaking parameters related to the masses and mixings of the superpartners, 
which are not measured. It is impossible to scan the full parameter space of the general 
MSSM for phcnomcnological purposes. On the other hand, all of the soft parameters are 
presumably generated by some underlying supersymmetry breaking and mediation mecha- 
nism. Usually, such a mechanism will provide relations between the MSSM soft parameters 
and reduce the number of parameters significantly. The properties make it possible, after 
the detection of low energy supersymmetry, to distinguish different supersymmetry breaking 
scenarios based on experimental data. Since the LSP appears naturally in the MSSM as 
a natural candidate for cold dark matter, it is interesting to study the prospects for dark 
matter observations in different SUSY breaking scenarios. 

Most of the scenarios we consider in this article fall into the general category of gravity 
mediated supersymmetry breaking [48]. Gravity is a natural candidate to mediate supersym- 
metry breaking because regardless of other details of the model, the gravitational interaction 
between the hidden and the observed sector exists. For some other mediation mechanism to 
dominate, one must suppress the contribution from gravity. For this reason, many bench- 
mark scenarios have been constructed and studied in the gravity mediation framework. The 
main alternative to gravity mediation is supersymmetry breaking mediated by gauge inter- 
actions [49]. In the gauge mediation scenario, we have approximately the following relation 
between the soft masses and the gravitino mass 

— ~ ^ < 1, (4) 

msoft «a Mpi 

where Ms is some typical supersymmetry breaking scale. Therefore, generically, we will 
have a very light gravitino as the LSP [50]. Since we focus our attention in the paper on 
neutralino dark matter, we will not consider gauge mediated scenarios in detail. 

Before discussing the detailed study, we briefly review the different scenarios we have con- 
sidered. Notice that the patterns of soft parameters resulting from different SUSY breaking 
scenarios are given at the input scale, which we take to be the grand uniflcation scale. The 
phenomenologically interesting low energy parameters are obtained through the running of 
the renormalization group equations. 

A. mSUGRA 

The general framework to study gravity mediated supersymmetry breaking is the N=l, 
D—A supergravity Lagrangian which is an effective Lagrangian after integrating out quantum 
gravity effects at the Planck scale. Without a specific model, general gravity mediation will 
not give rise to specific predictions and relations between different parameters. The most 
general supergravity Lagragian would contain all of the terms allowed by gauge symmetry. 
However, a minimal benchmark in this scenario was constructed [48] and well studied as an 
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example of the physics of gravity mediated supersymmetry breaking. The main ingredients 
of this scenario (at the supersymmetry breaking scale) are 

1. Flavor diagonal: 

The left-left and right-right blocks of sfermion mass matrices are diagonal and pro- 
portional to the unit matrix. The trilinear couplings are also universal. This type of 
sfermion mass matrix is primarily motivated by the phenomenological consideration 
of satisfying flavor constraints. 

2. Universahty: 

All of the diagonal soft masses are assumed to have the same value, itlq. The diagonal 
trilinears are also assumed to be the same, with value Aq. This is a assumption 
motivated largely by the simplicity of the parameter space. 

3. Universahty of the gaugino masses, M1/2: 

Gauge unification, in simple GUT scenarios, implies the unification of the gaugino 

masses. Therefore, the mSUGRA choice of gaugino mass universalilty is primarily 
motivated by gauge coupling unification. Notice however, the link between these two 
issues is not always as direct as in the simplest scenario. Hence, as with the requirement 
of gauge unification, we should still treat gaugino mass unification as an assumption. 

Under those assumptions, the mSUGRA model contains the following set of parameters 

mo, ^0, Ml/2, tan/3,Sign(//). (5) 

In most of the parameter space of the mSUGRA model, the LSP is Bino-like (the super- 
partner of the hypercharge gauge boson) . 

We note that despite its simphcity, mSUGRA only represents a very special corner of the 
possible model space of supersymmetry breaking. Most of the universality assumptions may 
not be present in other scenarios while still satisfying the motivations of gauge unification 
and respecting fiavor constraints. 

B. Anomaly Mediated Supersymmetry Breaking 

The authors of Ref. [51] pointed out an important source of supersymmetry breaking 
mediation in the supergravity Lagrangian due to the so-called super-conformal anomaly. 
It produces a elegant solution to the flavor problem due to the "UV-insensitivity" of the 
soft parameters in the sense that they can be expressed purely in terms of the low energy 
parameters such as the Yukawa couplings and gauge couplings. Reviewing the details of the 
origin of such a source is outside of the scope of this paper. Instead, we will focus on the 
characteristics of the soft spectrum as a result of anomaly mediation. The soft spectrum of 
the minimal anomaly mediation is [51,52] 

9a 
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= -^^3/2, (6) 

where y collectively denotes the Yukawa couplings. The /3-functions and anomalous dimen- 
sions, 7, are computed in the supersymmetric limit [51,52] and, therefore, are functions of 
the gauge couplings and superpotential parameters, mo is some addition parameter outside 
of the AMSB scenario, not to be confused with the universal scalar mass in the SUGRA 
scenario. It is added by hand mainly to give rise to a non-tachyonic stau mass. In the 
minimal AMSB scenario, the scalar masses squared at the GUT scale are each given by: 




^(3y ) ml/2 + D terms. (7) 



However, in non-minimal scenarios, the first term may vary. 

The most obvious feature of the anomaly mediation is the proportionality of the (3- 
f unction to the soft parameters. Since the origin of the mediation mechanism is the super- 
conformal anomaly, or the scale dependence of the parameters of the Lagrangian, it is not 
surprising to see this proportionality to the slope of the RGE running of the parameters. 
Notice some important features of the anomaly mediation spectrum: 

1. The main phenomenological implications of this scenario are that the gaugino masses 
have the ratios: 

Ml : M2 : Ms = 2.8 : 1 : 7.1 (8) 

leading to the LSP being a neutral wino [52,53], which is only slightly lighter than the 
charged wino (by a few hundred MeV). This leads to a long lived lightest chargino 
with a distinctive collider signature. Notice that in the AMSB scenario, the gluino 
is usually quite heavy due to its large input scale value. This will in turn lead to 
large squark masses, since typically m| ~ com|(0) -|- C3M|(0), where cq ~ 0(1) and 
C3 ~ 5 — 7. 

2. Unfortunately, the slepton masses squared turn out to be negative due to the fact 
that the pure anomaly mediation contribution to its value is negative. This is clearly 
unacceptable as it will lead to a charge breaking minima. This is also the reason 
behind the introduction of the parameter mo as it parameterizes additional physics 
which may solve this problem. However, introducing such a parameter undermines one 
of the main advantages of the AMSB scenario. Since there is no good reason to assume 
the additional physics is universal, this may reintroduce the flavor problem. Solutions 
to this problem have been proposed [54], for example by combining this mechanism 
with an anomalous U{1) gauge group [55]. 

AMSB is a very characteristic corner of the model space of the supersymmetry breaking. 
It could be important if the tree level gravity mediation (which is generically one-loop order 
larger than the AMSB piece) is suppressed. If AMSB is indeed the dominant contributor 
to the SUSY breaking, the suppression of other potential contributions (which are not sup- 
pressed in a generic supergravity Lagrangian) will be an important clue to the structure of 
the fundamental theory. 
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C. Other benchmarks 



Although mSUGRA and AMSB are two of the most widely studied benchmark scenarios, 
there are other benchmarks proposed in recent years with different motivations. Here we 
give a brief survey of a few of them which we will consider. 

1. Focus Point 

Focus point supersymmetry [56], in its simplest reahzation, is a special corner of the 
mSUGRA parameter space. In this corner, very large scalar masses are possible without 
violating naturalness constraints. The resulting low energy parameters have a very important 
feature. Namely, the soft masses squared of the Higgs boson, mjj^, very important in 
electroweak symmetry breaking, have pseudo fixed point behavior at the electroweak scale. 
They can start with a wide range of input values and run to a similar negative value at the 
low scale. This is interesting because it indicates that, in the focus point region, electroweak 
symmetry breaking does not require fine-tuning in the high energy input values. 

A typical feature of the focus point region is large scalar soft masses (usually ~ TeV). 
The LSP is usually a higgsino or a mixed higgsino-gaugino. This is quite different from 
the other cases where the LSP has a dominant gaugino content. The main reason for a 
larger higgsino content is again the larger input value of the soft scalar mass. The tree level 
electroweak symmetry breaking condition gives 

-ruz rriH^ - // . (9) 

In the typical mSUGRA scenarios, mf^^ is driven to some large negative value due to the 
RGE running. This requires a large value of fi to give the correct Z mass. However, in the 
focus point region, it is possible that the large input value of the scalar soft mass makes 
less negative. Hence, a smaller value of /i is possible, which leads to a larger higgsino 
content in the LSP. 

2. Michigan Benchmarks 

Recently, another set of benchmark models has been proposed [57]. The main motivation 
of these scenarios is to consider different patterns of gaugino masses (compared to the 
standard relation in mSUGRA, for example). By doing so, the gluino mass does not need 
to be quite so heavy (as implied by mSUGRA or AMSB relations). It is shown in Ref. [57] 
that the fine-tuning of the electroweak symmetry breaking in the MSSM is dominated by 
the gluino mass. Therefore, a different prediction for the gluino mass will generically offer 
hope of reducing fine-tuning. 

Another distinctive feature of this set of benchmark models is that each of them has 
a motivation in string theory, rather than making arbitrary simplifying assumptions about 
the pattern of the soft parameters at the input scale. 

Most of these proposed benchmarks fall in the category of gravity mediated supersym- 
metry breaking. The main feature of these constructions is the suppression of the tree level 



10 



dilaton contribution which give rise to mSUGRA-hke universal input values. Therefore, 
they are essentially different combinations of 1-loop contributions such as anomaly medi- 
ation, threshold corrections, and Green-Schwarz counter terms. There are also two cases 
in which the spectrum resembles the typical gauge mediated scenario but without a light 
gravitino. 

For more phenomenologically minded readers, the most important information is the 
spectrum of superpartners presented in tables 1 and 2 of Ref. [57]. 

3. Gaugino Mediation 

Gaugino mediated supersymmetry breaking, proposed in Ref. [58], represents another 
class of SUSY breaking mediation motivated by the brane-world scenario. It achieves the 
suppression of unwanted supersymmetry breaking effects, such as the flavor violating cou- 
plings, by the separation of the observable and hidden sectors via the separation of their 
respective branes. In the simplest realization, it begins with a two branc set up. Most of 
the matter flelds, in particular the chiral fermions and their superpartners, localize on one 
of those branes while the supersymmetry breaking sector localizes on the other. If this is 
the whole description, it is again similar to the anomaly mediation scenario. The new idea 
in this scenario is that the vector multiplet, in particular, the gauginos, are now allowed to 
propagate in the bulk. 

The pattern of soft masses in this scenario is straight forward to describe: 

1. The gauginos will have "direct" couplings to the supersymmetry breaking sector. 
Therefore, their soft masses are proportional to F/M, where F is supersymmetry 
breaking order parameter and M is some underlying fundamental scale which char- 
acterizes the coupling between gauginos and the hidden sector (since the coupling is 
usually of the form of a non-renormalizable term suppressed by M). With proper 
choices of F and M, the gaugino masses in this scenario can be chosen to be similar 
to any of the other supersymmetry breaking mediation scenarios. 

2. The soft masses in this scenario cannot come from the direct couplings due to the 
separation of the two branes. On the other hand, they can be generated from the 1-loop 
diagrams in which a gaugino is emitted, travels through the bulk to the supersymmetry 
breaking brane, gets the information of SUSY breaking and then returns to join the 
sfermion propagator again. Generically, the soft masses are then suppressed relative 
to the gaugino mass by a loop factor, ~ M|/(167r^). 

3. At low energy, soft masses and trilinear couplings generated through the RGE runnings 
are all proportional to the gaugino masses and are universal (since gauge coupling 
is flavor diagonal). Therefore, in this scenario, any flavor violating effect is again 
suppressed. 

Notice that the soft masses for sfermions are always positive in this model since they 
essentially start from zero at the high scale and get positive corrections from RGE running. 

In this paper, we have considered only a small sample of possible gaugino mediated 
models in which the features of mSUGRA are manifest while respecting the conditions of 
point 2 (m^-- M|/(167r2)). 
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V. CALCULATION OF OBSERVABLE QUANTITIES AND PARAMETER SCAN 



Assuming R-parity is conserved, a large portion of the parameter space in the scenarios 
mentioned above can provide a neutrahno which is a viable dark matter candidate. The 
present relic density for such a particle can be calculated by solving the Boltzman equa- 
tions. This calculation is more complicated, however, if the mass of another supersymmetric 
particle is only slightly greater than that of the LSP. In this situation, these particles may 
efficiently co-annihilate^, lowering the rehc abundance of a thermal rehc [59]. We have used 
the DARKSUSY package for our rehc density calculation which includes all coannihilations 
with charginos and neutralinos as well as all resonances and thresholds in the cross sec- 
tion calculation [60]. To calculate the supersymmetric particle spectrum, we have used the 
packages SUSPECT [61] and FeynHiggs [62]. 

The thermal rehc density of neutrahnos can be compared to the WMAP result of 
^CDM^^ = O.llSloiois [4]- It is possible, however, that the present density of dark mat- 
ter particles may have been produced by non-thermal mechanisms [52,63]. With this in 
mind, we treat the WMAP dark matter density measurement as an upper limit for the ther- 
mal relic abundance of any given WIMP. For models which predict a thermal relic density 
below the WMAP measurement, we assume that non-thermal mechanisms are sufficient to 
generate the measured dark matter density. 

There can be strong constraints on supersymmetric models imposed by measurements 
of the 6 — S7 branching fraction [64]. In anomaly or gaugino mediated scenarios, the 
supersymmetric flavor problem is solved naturally and this constraint is typically not vio- 
lated. In the general MSSM or mSUGRA, however, many of the models violate this bound. 
On the other hand, we notice that one should treat these constraints with caution since a 
small deviation from mSUGRA, such as a small off diagonal term in the squark mass ma- 
trix, could change significantly the prediction for h ^ without effecting the dark matter 
phenomenology significantly. 

Ref [61] also present a scenario, called the "phenomenological" MSSM, which covers a 
subset of the full MSSM parameter space (18 free parameters). It essentially relaxes the 
universality conditions in mSUGRA between the first two generations and the third, the 
two Higgs soft masses squared, and the three gaugino masses. All of the off diagonal entries 
in the sfcrmion mass matrices as well as the CP violating phases are still set to be zcro^. 
Although this "truncated'' parameterization of the MSSM is not specifically linked to any 
supersymmetry breaking mechanism and is not general enough to give a model independent 
conclusion, it is useful to illustrate what the dark matter signature may look like if the real 
supersymmetry breaking scenario is outside of the models we have considered. Therefore, 



Co-annihilations are annihilations between LSPs and other supersymmetric particles, such as 
heavier neutralinos, charginos, staus, etc. 

^Notice that CP violating phases could have important impact on the dark matter phenomenology. 
For example, they could significantly change the gaugino-quark-squark coupling and thereby have 
a significant impact both on annihilation and elastic scattering cross sections [65]. 
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we also include this more general scenario in our study. 

At the 2001 Snowmass meeting, a set of benchmark supersymmetry scenarios were pro- 
posed [66]. These benchmarks include representative models within mSUGRA, as well as 
gauge and anomaly mediated models. These selected cases help to illustrate the various phe- 
nomenological features associated with different areas of the large MSSM parameter space. 
Although in principle they are already included in the parameter scan of mSUGRA and 
AMSB, it is interesting to consider them separately as representative corners of the param- 
eter space. Therefore, we have studied the dark matter phenomenology for many of these 
selected benchmarks and included them in our figures. 

VI. RESULTS AND DISCUSSION 



A. Results For Direct Detection Experiments 

The rates and sensitivities for direct detection experiments can be reliably calculated 
with relatively small uncertainties which result from measurements of the local dark matter 
density. For this reason, results from direct detection experiments are among the most useful 
tests of supersymmetric dark matter. 

In figures 1 and 2, we show comparisons of the neutralino mass to the scalar (spin inde- 
pendent) neutralino-nucleon elastic scattering cross section for a variety of models. In figure 
1, the lightly shaded area corresponds to a the phenomenological MSSM search described in 
the previous section. The darker shaded region represents those models limited to minimal 
supergravity (mSUGRA). For this case, we have allowed mo, Tni/2 and Aq to vary between 
and 10,000 GeV. We selected tan/3 to have values in the range of 1 to 50 and have allowed 
to take on either positive or negative values. We acknowledge that 10,000 GeV is probably 
outside of the interesting range for mi/2 or Aq, but we choose to allow such values in an 
attempt to be inclusive. Our results would vary only slightly if we were to use an upper 
limit of a few TeV, for example. Also shown as dark circles are the points which we consider 
to be gaugino mediated supersymmetry breaking models. In such models, mo is less than 
47rmi/2- In this manifestation, gaugino mediation is a subset of mSUGRA. 

For each point shown, the relic density is below the maximum value allowed by the 
WMAP data (fi^/i^ < 0.129). Models that violated accelerator hmits are also not shown, 
except for b to S7 limits which may be violated. The experimental limits and sensitivities, 
shown as lines, from top to bottom (on the right) are CDMS-March 2002 (solid) [22], ZEPLIN 
1-final 2002 (dashed) [23], Edelweiss-2000+2002 (dots) [24], CRESST Il-projected hmit 
(solid) [26], CDMS-projected limit (dots) [27], Edelweiss Il-projected hmit (dashed) [28], 
ZEPLIN 4-projection (sohd) [29] and XENON-1 ton projected hmit (dashed) [30]. Also 
shown, as a solid contour, is the region which the DAMA collaboration claims evidence [25]. 
For a summary of present limits and projections for all direct dark matter experiments, see 
Ref. [67]. 

Figure 1 indicates that very little of mSUGRA space has been effectively probed by 
existing direct searches. However, future experiments such as GENIUS and XENON will 
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have the abihty to probe nearly all mSUGRA models with LSP masses below about 200 
GeV. 

In figure 2, again the lightly shaded area corresponds to the phenomenological MSSM. 
The darker shaded region shows those models limited to minimal Anomaly Mediated Super- 
symmetry Breaking (mAMSB) models. In such models, we allowed mo to vary between 
and 10,000 GeV and the gravitino mass, m3/2, to vary between 25,000 and 3,000,000 GeV. 
Again tan j3 was selected between 1 and 50 and yU is allowed to have either sign. Wc acknowl- 
edge that 10,000 GeV or 3,000,000 GeV are probably outside of the interesting range for mo 
or m3/2, respectively, but we choose to allow such values in an attempt to be inclusive. Our 
results would vary only slightly if we were to use a smaller upper limit. 

We assumed universality of the scalar masses in the minimal AMSB model. Also shown 
in figure 2 are black X's which represent randomly selected non-minimal AMSB models. In 
these models, in addition to the parameters varied in the minimal AMSB, the scalar masses 
squared were allowed to vary from at the GUT scale. Each scalar mass squared was set 
to mQ multiplied by a factor randomly selected between 5 and -5. 

Figure 2 shows that even with future experiments, AMSB scenarios are very difficult 
to test with direct detection methods. In non-minimal models, these prospects are slightly 
better. A primary reason direct detection is difficult in AMSM models is that the value of 
M3 is considerably higher than M2 or Mi, and, therefore, squark masses tend to be much 
larger than the LSP mass. This leads to suppression in the s-channel diagram for neutralino- 
quark scattering. Additionally, AMSB typically predicts a gaugino-like LSP, less favorable 
for direct detection than a mixed gaugino-higgsino or pure higgsino. 

Figure 3, illustrating the effect of scalars in direct detection, shows the mass of the 
lightest neutralino versus the universal scalar mass, mo, in the mSUGRA scenario. The 
darkest points are those currently excluded by either CDMS-March 2002, ZEPLIN 1-final 
2002 or Edclwciss-2000+2002. The intermediate points arc testable by planned experiments 
including CRESST II, CDMS (Soudan), Edelweiss II or ZEPLIN 4. The hghtest points are 
beyond the sensitivity of ZEPLIN 4. 

Figure 4 shows the sensitivities to direct experiments for a variety of Snowmass and 
Michigan benchmark scenarios. Snowmass slopes la, 3 and 9 and shown as lines from left 
to right, respectively. Snowmass point lb, and samples of slopes 2, 4 and 5, are shown as 
b, 2, 4 and 5 in the figure. For the Michigan benchmark scenarios, models A through G are 
shown. 

Along the length of each slope shown (in the Snowmass plot), the relic density is below 
the maximum value allowed by the WMAP data {Vt^h"^ < 0.129). However, some of the 
points (rather than slopes) shown produce a larger relic density. 

B. Results For Gamma- Ray Experiments 

In the case of a cuspy or spiked galactic halo, gamma-ray astronomy can provide a 
strong test of supersymmetric dark matter. For models with flat galactic cores, however, 
such observations are nearly impossible in planned experiments. 

In figure 5, the rate of continuum gamma-rays above 1 GeV per square meter, per year 
of exposure, from the galactic center is shown verses neutralino mass. A smooth NEW halo 
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profile is assumed. If a Moore, et. al. model were considered, each point would produce 
approximately a factor of 10^ more events; considerably more if a spiked halo model were 
used. As in figures 1 and 2, the lightly shaded region represents the phenomenological 
MSSM, the darker region corresponds to mSUGRA models, and the darker fine represents 
the AMSB models (minimal and non- minimal) . Also shown are limits from EGRET (solid) 
and the predicted sensitivity for the future experiment GLAST (dashed) [14]. Note that 
EGRET is already sensitive to some models, and GLAST will be sensitive to many models, 
especially for those models with a relatively light neutralino. In AMSB scenarios, the mass 
parameters Mi and M2 are proportional to their respective couplings, so the neutralino mass 
and annihilation cross section become correlated, leading to the single line across the figure. 

Figure 6, ilhistrating the effect of scalars in gamma-ray production, shows the mass 
of the lightest neutralino versus the universal scalar mass, mo, in the mSUGRA scenario. 
The darkest points are those currently excluded by EGRET assuming a smooth NEW halo 
profile. The intermediate points are testable by the GLAST experiment. The lightest points 
are beyond the sensitivity of GLAST. 

Again, we also show results for Snowmass and Michigan benchmark scenarios. Figure 7 
shows the sensitivities to gamma-ray experiments for these models. The models shown are 
the same as those shown for the direct experiments. 

In figures 8 and 9, the rate of gamma-rays from the line producing processes, x°X° ~^ 11 
and XX. ~^ 1^, are shown [15]. The sensitivity of GLAST is shown (dashed) for comparison. 
Again a smooth NEW halo profile is assumed. AMSB scenarios are the best prospect for 
line detection, although if a more cuspy model, such as Moore, et. al., were realized, even 
some mSUGRA models may be observable. 

C. Results For Neutrino Experiments 

Neutrinos from WIMP annihilation in the Sun (or Earth), like direct detection, provide a 
method of dark matter detection with rates and sensitivities which can be reliably calculated. 

In figures 10 and 11, we show the rates in a kilometer scale neutrino telescope, such as 
IceCube, per year from neutralino annihilation in the Sun. The models are represented by 
the various types of shading as in figures 1 and 2. The sensitivity projected for the IceCube 
experiment (dashed) is also shown [68] . A 50 GeV muon energy threshold has been imposed. 

Some models in the phenomenological MSSM with a substantial higgsino fraction can 
have very large scattering cross sections and, therefore, provide observable neutrino rates. 
On the other hand, in either mSUGRA or AMSB scenarios, few models can be probed with 
kilometer scale neutrino telescopes. An important exception is the foucs point region of the 
mSUGRA parameter space, in which the LSP also has a large higgsino component. These 
rates do not depend strongly on the choice of halo profile as is the case for observations of 
the galactic center. 

Again, we also show results for Snowmass and Michigan benchmark scenarios. Figure 12 
shows the sensitivities of neutrino experiments for these models. The models shown are the 
same as for the direct and gamma-ray experiments. 
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D. Results For Positron and Anti-Proton Experiments 



Figure 13 compares the positron boost factor required to explain the HEAT data to the 
maximum anti-proton boost factor which does not violate the BESS bound. Results are 
shown for a variety of models and SUSY breaking scenarios. The models are represented 
by the various types of shading as in figure 5. It is somewhat difficuh to interpret these 
results, however, due to large uncertainties in the galactic halo model. First, the question of 
how large the positron boost factor could be is a subject of some debate. It is known that 
on the very large scales of galaxies and clusters of galaxies, such enhancement factors can 
exceed ~ 10^ according to simulations [69]. Perhaps boost factors significantly above this 
scale could be considered unrealistic. 

A second uncertainty is the degree to which the positron and anti-proton boost factors 
could differ. A ratio of 1 to 10 between the two factors is certainly reasonable, but it is 
difficult to determine how much larger this ratio could be. In figure 13, lines are shown 
representing the cases of equal positron and anti-proton boost factors and a positron boost 
factor ten times the anti-proton boost factor. Perhaps models to the upper-left of this second 
line are less likely to be capable of explaining the observed positron excess. 

Note that even for a smooth halo profile (boost factors of one) , a fraction of the super- 
symmetry models considered produced an anti-proton fiux larger than would be consistent 
with the BESS data. 

An LSP with lesser Bino content, such as in the AMSB (wino) or focus point (higgsino) 
scenarios, will have a greater chance of producing the observed positron excess. Such an 
LSP may annihilate dominantly into W bosons (if kinematically allowed) which makes the 
positron spectrum considerably easier to fit. Also, not that our positron analysis is not a fit 
to the observed complete positron spectrum, but rather only a fit to a single point (near the 
observed peak), disregarding the spectral shape. We also assume the background positron 
flux is known [70]. A more careful analysis should be done using a fit to the entire 
observed spectrum, with the normalization of the background positron flux treated as a free 
parameter. The result of such an analysis [18] will be that there are supersymmetric models 
capable of producing the excess without a prohibitively large boost factor. Our analysis 
presented here is qualitatively right in the sense that it compares the relative potential for 
producing a positron excess between different models, which is what we are focusing on here. 

E. Correlations Between Different Experimental Searches and Constraints 

In some cases, the neutralino characteristics which determine the observability of one 
detection method also play an important role for another method. It is useful to identify 
these features in order to discuss the prospects for neutralino dark matter detection over 
multiple experimental methods. 

Prospects for both direct detection experiments and indirect detection experiments which 
rely on WIMP scattering on the Sun or Earth depend on the elastic scattering cross section 
of WIMPs with nucleons. For this reason, the fiux of neutrinos from the Sun from neutralino 
annihilation is related to the event rate in direct detection experiments. Figure 14 shows the 
rate of muons (from charged current neutrino interactions) from neutralino annihilation in 



16 



the Sun for the phenomenological MSSM. Black points represent models which have allready 
been excluded by current experiments such as CDMS-March 2002 [22], ZEPLIN 1-final 2002 
[23] and Edelweiss-2000+2002 [24]. Darker points represent those models which can be 
tested by planned experiments with a sensitivity near that of the ZEPLIN 4-projection [29]. 
Lighter points fall below this sensitivity. 

There is a clear correlation between the sensitivities of these two types of experiments. 
For models with a neutralino mass substantially above the energy threshold of neutrino 
telescopes, the sensitivity of an experiment such as IceCube is very similar to that of next 
generation direct detection experiments. Given this observation, these two methods of de- 
tection will likely be capable of confirming an observation made by the other technique. 
Note that this relationship is not present for all SUSY breaking scenarios, however. 

Several methods of indirect detection depend strongly on the neutralino-neutralino an- 
nihilation cross section. For this reason, event rates for these experiments arc correlated. 
For example, in figure 15, we show the rate of gamma-rays above 1 GeV, per square me- 
ter, per year of exposure, from the galactic center, verses neutralino mass, as in figure 5. 
Particularly high neutralino annihilation cross sections will not only result in a high rate 
of gamma-rays, but also a high rate of other annihilation products, including anti-protons 
from the galactic halo. Models in the lightly shaded region do not overproduce anti-protons 
in the case of a smooth halo profile, as measured with the BESS experiment [47]. Models in 
the darker region, however, require some degree of fine tuning in the halo model, i.e. placing 
clumps at large distances, to accomodate this limit. The consequence of this comparison 
is the conclusion that for those models in which EGRET is sensitive, they are also disfa- 
vored by anti-proton measurements. If, however, a more cuspy halo model is considered, the 
rate from annihilations in the galactic center region will increase, although the anti-proton 
and positron rates, which depend primarily on the local halo characteristics, will not be 
significantly affected. 

If constraints from the measurements of the b sj branching fraction are considered, 
many models in the mSUGRA scenario can be strongly disfavored. In figures 16 and 17, we 
show the impact of the 6 — > 57 constraint on the prospects for direct detection. In figure 16, 
all sample points shown correspond to the phenomenological MSSM, as described earher. 
Figure 17 shows sample points in the mSUGRA scenario. In each figure, lighter points do 
not violate the b 57 constraint. Darker points do violate this constraint. 

Note that those models most easily tested by direct detection experiments are more 
likely to violate the 6 — > 57 constraint. This is especially apparent in the case of mSUGRA 
where, except for small pockets, the only models which do not violate b ^ sj are beyond 
the sensitivities of all present and next generation direct experiments. One pocket in which 
this constraint is not violated, near 60-100 GeV and 10^^ pb, are those mSUGRA points 
which are gaugino mediated. Such features make scenarios such as gaugino mediation and 
AMSB, which do not generally violate this constraint, appear very attractive. 

F. Comparison of Diflferent Supersymmetry Breaking Scenarios 

From our numerical results, we have seen that different supersymmetry breaking sce- 
narios generically give rise to different dark matter signatures since they predict different 
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TABLE L A comparison of the detection potential and consistency of different SUSY break- 
ing models. For direct searches, continuum gamma-rays from the galactic center, line emission 
gamma-rays from the galatic center and neutrinos from the sun, each breaking scenario is de- 
scribed as "Good", "Fair" or "Difficult". These assessments reflect the likelihood of detection 
with a given method. For the positron column, "Likely" , "Less Likely" and "Unlikely" reflect the 
likelihood of neutralino annihilation in the galactic halo being capable of producing the observed 
positron excess. Finally, the anti-proton column evaluates whether the predicted anti-proton flux is 
compatiable with the BESS measurement, or if a fine tuned halo distribution is needed to reconcile 
such a model. 
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LSP properties as well as different spectra for the other superpartners. We qualitatively 
summarize these differences in Table I. 

In Table I, we hst the mSUGRA and AMSB models first, as they are the best studied 
scenarios. We then grouped the michigan benchmarks into four groups based on the identity 
of the LSP. The gaugino mediation scenario and the Snowmass benchmark points are also 
listed. AlthoTigh these points are allready covered in the scan of mSUGRA and AMSB, we 
find it illustrative to study these benchmarks seperately as well. 

We then assign qualitative measures such as "good", "difficult", etc., to each super- 
symmetry breaking scenario regarding their potential to produce observable signals in a 
particular experiment. For the positrons column, we do not attempt to assess the liklihood 
of observing such a signal, but rather estimate the plausibility that the observed positron 
excess could be generated from WIMP annihilations for a given SUSY model. For the anti- 
protons column, we assess whether the SUSY model requires fine tuning of the halo profile 
to remain consistent with the BESS data. 

Notice that mSUGRA is not a point but rather a scenario with a sizable parameter 
space. Therefore, the qualitative comments assigned to mSUGRA should be understood as 
a general assessment applied to much of the parameter space. Different corners of mSUGRA 
parameter space could have quite different dark matter phenomenology. The same is true 
for our assessment of AMSB. Notice that these comments are necessarily qualitative and 
are not objective measures. For a more precise assessment, see the figures we present in this 
paper rather than the table. 

A main lesson we can take away from Table I is that even if a single experiment cannot 
provide a definitive answer to distinguish between different supersymmetry breaking sce- 
narios, a combination of the results of all the experiments could be a powerful tool in this 
respect. For example, several scenarios in this table could produce a sizeable signal in future 
direct search experiments. However, a observation (or non-observation) from the neutrino 
signal from the Sun/Earth could favor (or disfavor) the focus point region of the parameter 
space. Another example is the combination of the positron signal with direct detection. 
As we can see from the table, AMSB, Mich D,E, and the focus point region can give rise 
to a sizeable positron flux. However, their direct detection rates differ quite significantly. 
Notice that most of the models we study wiH at least give rise to a potentially observable 
signal in some experiment, and it will be possible to evaluate these scenarios based on the 
experimental results. 

VII. CONCLUSIONS AND OUTLOOK 

The large number of experiments designed to search for supersymmetric dark matter par- 
ticles are beginning to enter a very exciting phase. Direct, indirect and collider experiments 
are each capable of a great deal of progress in the near future. 

1. Direct Searches: 

Direct searches are expected to improve their sensitivity by more than a factor of 10^ 
with planned experiments such as ZEPLIN-4 [29], GENIUS [31], or a 1 ton version of 
XENON [30]. With such improvements, many scenarios, including much of mSUGRA, 
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will be probed by these experiments. AMSB scenarios are more difficult to test by 
direct experiments, although some models should be accessible. Gaugino mediated 
models may be testable in these planned experiments. 

2. Neutrinos Prom Solar Capture: 

Indirect experiments which involve the capture of neutralinos in the center of the Sun 
(or Earth) depend principally on the neutralino-nucleon elastic scattering cross section 
for capture and, therefore, have a tendency to be sensitive to the same models as direct 
experiments. In fact, we have shown that for moderate or heavy neutralinos (~ 150 
GeV or heavier), IceCube and ZEPLIN-4 have very similar sensitivities. The focus 
point region of mSUGRA is, perhaps, the best prospect for this experimental method. 

3. Gamma-Rays Prom The Galactic Center: 

Rates for gamma-ray (as well as positron and anti-proton) experiments depend pri- 
marily on the neutralino annihilation cross section. Por these methods, the galactic 
halo profile can be of enormous importance. Particularly, the sensitivity of gamma-ray 
experiments, such as EGRET or GLAST, depend critically on the halo profile. Cuspy 
halo profiles often provide observable gamma-ray rates in such experiments, although 
profiles with lower density cores generally do not. Continuum emission of gamma-rays 
would be most easily observed in AMSB scenarios, although the prospects are still 
quite good for mSUGRA models. Gaugino mediated models arc somewhat more dif- 
ficult. Por line emission, however, gaugino mediated models provide better rates than 
other mSUGRA scenarios. 

4. Positrons and Anti-Protons: 

Cosmic positron and anti-proton experiments do not depend strongly on the size of the 
cusp, but rather on the smoothness or dumpiness of dark matter in the local region, 
as well as the location of any nearby clumps. As was the case for gamma-ray rates, 
positron and anti-proton rates are largest in AMSB scenarios. 

5. Plavor Constraints: 

Plavor constraints, b 57 in particular, can reveal important information relevant to 
neutralino dark matter phenomenology. In breaking scenarios without a solution to the 
flavor problem, SUGRA for example, many models violate this constraint, especially 
those which are most easily probed by direct experiments, or indirect experiments 
which use WIMP capture in the Sun. Scenarios which include a solution to the flavor 
problem, such as AMSB, do not often have such features. 

6. Correlations and Comparisons: 

By combining the results of multiple experimental techniques, characteristics of a 
particular supersymmetry breaking mechanism may be observed. Por example, rates 
which depend primarily on the LSP annihilation cross section could be combined with 
rates which depend on the LSP-nucleon elastic scattering cross section. We find that 
by combining all the experimental data which can be obtained in the foreseeable future, 
it may be possible to actually make distinctions between a variety of SUSY breaking 
scenarios. 
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7. Other Future Prospects: 

With the current operation of the Tevatron and future running of the Large Hadron 
CoUider (LHC), it is very hkely that supersymmetry will be discovered. With measure- 
ments of the supersymmetric spectrum (or portions of it), the problem of predicting 
event rates for direct and indirect dark matter experiments will likely become much 
more simple. Additionally, astronomical measurements of the galactic halo profile will 
shed a great deal of light on the prospects for indirect detection. It appears that a 
great deal of information is beginning to converge and an answer to these questions 
may not be far away. 
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FIG. 1. The scalar (spin-independent) neutralino-nucleon elastic scattering cross section verses 
neutralino mass. Each point corresponds to a theoretical prediction for a specific SUSY model. 
The light shaded area corresponds to a general or phenomenological MSSM model. The darker 
colored region are those models limited to the constrained MSSM (mSUGRA). Also shown as dark 
circles are the points which wc describe as gaugino mediated supersymmetry breaking models. 
See the text for more details. For each point shown, the relic density is below the maximum 
value allowed by the WMAP data {Cl^h^ < 0.129). Models that violate accelerator limits are 
not shown, except for b to 57 limits which will be discussed later. The experimental limits and 
sensitivities, shown as lines, from top to bottom (on the right) arc CDMS-March 2002 (solid) [22], 
ZEPLIN 1-final 2002 (dashed) [23], Edelweiss-2000+2002 (dots) [24], CRESST Il-projected hmit 
(solid) [26], CDMS-projected limit (dashed) [27], Edelweiss Il-projected hmit (dots) [28], ZEPLIN 
4-projection (solid) [29] and XENON-1 ton, projected limit (dashed) [30]. Also shown, as a solid 
contour, is the region which the DAMA collaboration claims discovery [25]. For a summary of 
present limits and projections for all direct dark matter experiments, see Ref. [67]. 
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FIG. 2. The scalar (spin-independent) neutralino-nucleon elastic scattering cross section verses 
neutralino mass. Each point corresponds to a theoretical prediction for a specific SUSY model. 
The light shaded area corresponds to a general or phenomenological MSSM model. The darker 
shaded region shows those models limited to minimal Anomaly Mediated Supersymmetry Breaking 
(mAMSB) models. Also shown are black X's which represent non-minimal AMSB models. See 
the text for more details. For each point shown, the relic density is below the maximum value 
allowed by the WMAP data {il.^h'^ < 0.129). Models that violate accelerator limits are also 
not shown, except for b to 57 limits which will be discussed later. The experimental limits and 
sensitivities, shown as lines, from top to bottom (on the right) are CDMS-March 2002 (solid) [22], 
ZEPLIN 1-final 2002 (dashed) [23], Edclweiss-2000-F2002 (dots) [24], CRESST Il-projcctcd limit 
(solid) [26], CDMS-projected limit (dashed) [27], Edelweiss Il-projected hmit (dots) [28], ZEPLIN 
4-projection (solid) [29] and XENON-1 ton, projected limit (dashed) [30]. Also shown, as a solid 
contour, is the region which the DAMA collaboration claims discovery [25]. For a summary of 
present limits and projections for all direct dark matter experiments, see Ref. [67]. 
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FIG. 3. The mass of the Ughtest neutraUno versus the universal scalar mass, mo, in the 
mSUGRA scenario. The darkest points are those currently excluded by one or more of 
CDMS-March 2002, ZEPLIN 1-final 2002 or Edelweiss-2000+2002. The intermediate points are 

testable by planned experiments including CRESST II, CDMS (Soudan), Edelweiss II or ZEPLIN 
4. The lightest points are beyond the sensitivity of ZEPLIN 4. The models shown meet the 
requirements of the previous figures. 
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FIG. 4. The scalar (spin-independent) neutralino-nucleon elastic scattering cross section verses 
neutralino mass. Each point corresponds to a theoretical prediction for a specific Snowmass (left 
frame) or Michigan (right frame) benchmark model. In the left frame, Snowmass slopes la, 3 and 9 
are shown as lines from left to right, respectively. Snowmass points lb, and samples of 2, 4 and 5 are 
shown as B, 2, 4 and 5 in the figure. In the right frame, Michigan benchmark models A through G 
are shown. See the text for more details. Along the length of each Snowmass slope shown, the relic 
density is below the maximum value allowed by the WMAP data (O^/i^ < 0.129). However, some 
of the points (rather than slopes) shown produce a larger relic density. The experimental limits and 
sensitivities, shown as lines, from top to bottom (on the right) are CDMS-March 2002 (solid) [22], 
ZEPLIN 1-final 2002 (dashed) [23], Edelweiss-2000+2002 (dots) [24], CRESST Il-projected hmit 
(solid) [26], CDMS-projected limit (dashed) [27], Edelweiss Il-projected hmit (dots) [28], ZEPLIN 
4-projection (solid) [29] and XENON-1 ton, projected limit (dashed) [30]. Also shown, as a solid 
contour, is the region which the DAMA collaboration claims discovery [25]. For a summary of 
present limits and projections for all direct dark matter experiments, see Ref. [67]. 
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FIG. 5. The rate of continuum gamma-rays above 1 GeV per square meter, per year of exposure, 
from the galactic center, verses neutrahno mass. A smooth NFW halo profile is assumed. If a 
Moore, et. al. model were considered, each point would produce approximately a factor of 10^ more 
events. As in figures 1 and 2, the lightly shaded region is for the general or phenomenological MSSM, 
the darker region corresponds to mSUGRA models, the shaded line represents AMSB models 
(minimal and non-minimal) and black circles indicate gaugino mediated models. For each point 
shown, the relic density is below the maximum value allowed by the WMAP data (n^h? < 0.129). 
Models that violate accelerator limits are also not shown, except for b to 57 limits which will be 
discussed later. See the text for more details. Also shown are limits from EGRET (solid) and the 
predicted sensitivity for the future experiment GLAST (dashed) [14]. Note that EGRET is already 
sensitive to some models, and GLAST will be sentivity to many models, especially for those models 
with a light neutralino. In AMSB scenarios, the mass parameters Mi and M2 are proportional to 
their respective couplings, so the neutralino mass and annihilation cross section become correlated, 
leading to the single line across the figure. 
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FIG. 6. The mass of the hghtest neutrahno versus the universal scalar mass, mo, in the 
mSUGRA scenario. The darkest points are those currently excluded by EGRET assuming a 
smooth NFW halo profile. The intermediate points are testable by the GLAST experiment. The 
lightest points are beyond the sensitivity of GLAST. The models shown meet the requirements of 
the previous figure. 
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FIG. 7. The rate of continuum gamma-rays above 1 GeV per square meter, per year of exposure, 
from the galactic center, verses neutrahno mass. A smooth NFW halo profile is assumed. If a 
Moore, et. al. model were considered, each point would produce approximately a factor of 10^ 
more events. In the left frame, Snowmass slopes la, 3 and 9 are shown as lines from left to 
right, respectively. Snowmass points lb, and samples of 2, 4 and 5 are shown as B, 2, 4 and 5 
in the figure. In the right frame, each point corresponds to a theoretical prediction for a specific 
Michigan benchmark model. Models A through G are shown. See the text for more details. Along 
the length of each slope shown, the relic density is below the maximum value allowed by the WMAP 
data {n^h^ < 0.129). However, some of the points (rather than slopes) shown produce a larger 
relic density. Also shown are limits from EGRET (solid) and predicted sensitivity for the future 
experiment GLAST (dashed) [14]. 
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FIG. 8. The rate of gamma-rays from the process X^X^ ~^ 77 P^r square meter, per year 
of exposure, from the galactic center, verses neutrahno mass. A smooth NFW halo profile is 
assumed. If a Moore, et. al. model were considered, each point would produce approximately 
a factor of 10^ more events. As in figures 1 and 2, the lightly shaded region is for the general 
or phenomcnological MSSM, the darker region corresponds to mSUGRA models, the shaded line 
represents AMSB models (minimal and non-minimal) and black circles indicate gaugino mediated 
models. For each point shown, the relic density is below the maximum value allowed by the 
WMAP data {Q.^h? < 0.129). Models that violate accelerator limits are also not shown, except 
for b to S7 limits which will be discussed later. See the text for more details. Also shown is the 
predicted sensitivity for the future experiment GLAST (dashed) [15]. In AMSB scenarios, the mass 
parameters Mi and M2 are proportional to their respective couplings, so the neutralino mass and 
annihilation cross section become correlated, leading to the single line across the figure. 
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FIG. 9. The rate of gamma-rays from the process X^X^ ~^ 7-^ P^r square meter, per year 
of exposure, from the galactic center, verses neutrahno mass. A smooth NFW halo profile is 
assumed. If a Moore, et. al. model were considered, each point would produce approximately 
a factor of 10^ more events. As in figures 1 and 2, the lightly shaded region is for the general 
or phenomcnological MSSM, the darker region corresponds to mSUGRA models, the shaded line 
represents AMSB models (minimal and non-minimal) and black circles indicate gaugino mediated 
models. For each point shown, the relic density is below the maximum value allowed by the 
WMAP data {Q,~^h? < 0.129). Models that violate accelerator limits are also not shown, except 
for b to S7 limits which will be discussed later. See the text for more details. Also shown is 
predicted sensitivity for the future experiment GLAST (dashed) [15]. In AMSB scenarios, the 
mass parameters Mi and M2 are proportional to their respective couplings, so the neutralino mass 
and annihilation cross section become correlated, leading to the single line across the figure. 
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FIG. 10. The rate of muons (from charged current muon neutrino interactions) above 50 GcV 
per square kilometer, per year, from neutrahno annihilation in the sun, verses neutralino mass. 
As in figure 1, the lightly shaded region represents the general or phenomenological MSSM, the 
darker region corresponds to mSUGRA models and black circles indicate gaugino mediated models. 
For each point shown, the relic density is below the maximum value allowed by the WMAP data 
{Q.^h? < 0.129). Models that violate accelerator limits are also not shown, except for b to S7 limits 
which will be discussed later. See the text for more details. Also shown is the predicted sensitivity 
for the next generation neutrino telescope IceCube (dashed) [68]. 
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FIG. 11. The rate of muons (from charged current muon neutrino interactions) above 50 GeV 
per square kilometer, per year, from neutrahno annihilation in the sun, verses neutralino mass. 
As in figure 2, the lightly shaded region is for the general or phenomenological MSSM, the darker 
shaded region show those models limited to minimal Anomaly Mediated Supersymmetry Breaking 
(mAMSB) models. Also shown are black X's which represent non-minimal AMSB models. For 
each point shown, the relic density is below the maximum value allowed by the WMAP data 
{Q.^hP' < 0.129). Models that violate accelerator limits are also not shown, except for b to 57 limits 
which will be discussed later. See the text for more details. Also shown is the predicted sensitivity 
for the next generation neutrino telescope IceCube (dashed) [68]. 
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FIG. 12. The rate of muons (from charged current muon neutrino interactions) above 50 GeV 
per square kilometer, per year, from neutrahno annihilation in the sun, verses neutralino mass. In 
the left frame, Snowmass slopes la, 3 and 9 and shown as lines from left to right, respectively. 
Snowmass points lb, and samples of 2, 4 and 5 are shown as B, 2, 4 and 5 in the figure. In the 
right frame, each point corresponds to a theoretical prediction for a specific Michigan benchmark 
model. Models A through G are shown. See the text for more details. Along the length of each 
Snowmass slope shown, the relic density is below the maximum value allowed by the WMAP data 
{n^h^ < 0.129). However, some of the points (rather than slopes) shown produce a larger relic 
density. Also shown is the predicted sensitivity for the next generation neutrino telescope IceCube 
(dashed) [68]. 
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Anti— Proton Boost Factor 

FIG. 13. The positron boost factor required to explain the HEAT positron excess [45,46] versus 
the maximum anti-proton boost factor consistent with the BESS anti-proton cosmic ray measure- 
ment [47]. As in figures 1 and 2, the Ughtly shaded region is for the general or phenomenolog- 
ical MSSM, the darker region corresponds to mSUGRA models, the shaded line corresponds to 
AMSB models (minimal and non-minimal) and black circles indicate gaugino mediated models. 
For each point shown, the relic density is below the maximum value allowed by the WMAP data 
{riy^hP' < 0.129). Models that violate accelerator limits are also not shown, except for b to limits 
which will be discussed later. Sec the text for more details. Solid lines correspond to the case 
when the positron boost factor is equal to the maximum anti-proton boost factor and for the case 
when the positron boost factor ten times larger than the maximum anti-proton boost factor. See 
the text for further discussion. 
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FIG. 14. The rate of muons (from charged current muon neutrino interactions) above 50 GeV 
per square kilometer, per year, from neutrahno annihilation in the sun, verses neutralino mass. All 
models show were randomly selected from the general or phenomenological MSSM, as described in 
the text. Black X's represent models which have allready been excluded by current experiments 
such as CDMS-March 2002 [22], ZEPLIN 1-final 2002 [23] and Edclwciss-2000+2002 [24]. Darker 
points represent those models which can be tested by planned experiments with a sensitivity near 
that of the ZEPLIN 4-projection [29]. Lighter points fall below this sensitivity. Also shown is 
the predicted sensitivity for the next generation neutrino telescope IceCube (dashed) [68]. For 
each point shown, the relic density is below the maximum value allowed by the WMAP data 
(r^^/i^ < 0.129). Models that violate accelerator limits are also not shown, except for b to 57 limits 
which will be discussed later. See the text for more details. 
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FIG. 15. The rate of continuum gamma-rays above 1 GeV per square meter, per year of ex- 
posure, from the galactic center, verses neutrahno mass. A smooth NFW halo profile is assumed. 
If a Moore, et. al. model were considered, each point would produce approximately a factor of 
10^ more events. Models in the lightly shaded region do not overproduce anti-protons in the case 
of a smooth halo profile, as measured with the BESS experiment [47]. Models in the darker re- 
gion, however, require some degree of fine tuning in the halo model, i.e. placing clumps at large 
distances, to accomodate this limit. Each model shown corresponds to a point in the general or 
phenomenological MSSM, as described in the text. For each point shown, the relic density is below 
the maximum value allowed by the WMAP data (O^^/i^ < 0.129). Models that violate accelerator 
limits arc also not shown, except for b to 57 limits. Sec the text for more details. Also shown are 
limits from EGRET (solid) and predicted sensitivity for the future experiment GLAST (dashed) 
[14]. Note that those models in which EGRET is sensitive, anti-proton measurements appear to 
disfavor. 
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FIG. 16. The impact of the b ^ s'j constraint [64] on direct detection for models in the general 
or phenomenological MSSM. Lighter points are models which do not violate this constraint. Darker 
points do violate this constraint. Also shown are various experimental limits and sensitivities as 
in figures 1 and 2. See the text for more details. 
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FIG. 17. The impact of the b ^ S'j constraint [64] on direct detection for models in the 
mSUGRA scenario. Lighter points are models which do not violate this constraint. Darker points 
do violate this constraint. Note that the small cluster of lighter points near 60-100 GeV and 10~^ 
pb are those mSUGRA points which we consider gaugino mediated (m^ ~ M|/(167r^)). Also 
shown are various experimental limits and sensitivities as in figures 1 and 2. See the text for more 
details. 
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